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The local-density approximation combined with the dynamical mean-field theory (LDA-I-DMFT) is 
applied to the paramagnetic phase of light actinide dioxides: UO 2 , Np02, and PUO 2 . The calculated 
band gaps and the valence-band electronic structure are in a very good agreement with the optical 
absorption experiments as well as with the photoemission spectra. The hybridization of the actinide 
5f shell with the 2p states of oxygen is found to be relatively large, it increases the filling of the 
5f orbitals from the nominal ionic configurations with two, three, and four electrons to nearly 
half-integer values 2.5, 3.4 and 4.4. The large hybridization leaves an imprint also on the core-level 
photoemission spectra in the form of satellite peaks. It is demonstrated that these satellites are 
accurately reproduced by the LDA-I-DMFT calculations. 

PACS numbers: 71.20.Gj, 79.60.-i, 71.15.-m, 71.27.-l-a 


I. INTRODUCTION 

Actinide dioxides are correlations-driven insulatorJ^^ 
that display a variety of complex ordered phases at low 
temperatures.!^ They crystallize in the CaF 2 structure 
(space group Fm3m), with eight-coordinated actinide 
atoms, and four-coordinated oxygen atoms. Due to the 
interplay of electron correlations, spin-orbital coupling 
and crystal field effects, the theoretical modeling of the 
electronic structure of these oxides presents numerous 
challenges. 

The conventional Kohn-Sham density-functional theory 
(DFT) in the local-density (LDA) and generalized gradient 
approximations fails to explain the insulating character 
of the oxides.!^ The band-gap problem was addressed a 
number of times using orbital-dependent functionals such 
as the self-interaction corrected LDA,!^ LDA -bU,!^ or the 
hybrid exchange-correlation functionals.®^ All of these 
calculations lead to insulators but the opening of the gap 
is intimately linked with the appearance of a long-range 
magnetic order. That is not satisfactory since the oxides 
retain the gap also in the high-temperature paramagnetic 
phase. Moreover, the orbital-dependent functionals predict 
a large magnetic moment at the plutonium atoms in PUO 2 
in disagreement with experiments.® These issues appear 
to be a rather general shortcoming of static mean-held 
approximations that build on a single determinant of 
Kohn-Sham orbitals. 

The dynamical mean-field theory (DMFT) is able to de¬ 
scribe correlated nonmagnetic insulatorsThis method, 
in combination with the density-functional theory, wa s suc- 
cessfully applied to selected actinide dioxides recently, 
and it indeed yields an insulating electronic structure with¬ 
out any long-range order and without any local magnetic 
moments in Pu02. In the present paper, we follow up on 
our previous study of the plutonium dioxide,where we 
employed a crystal-field potential deduced from experi¬ 


ments and assumed a simplified spherically symmetric 
hybridization of the plutonium 5f shell with the surround¬ 
ing electronic states. Here we relax these simplifications 
and determine the quantities entirely from first principles. 
We investigate also the paramagnetic phases of UO 2 and 
Np 02 in order to visualize trends in the behavior of the 
computed properties when the filling of the actinide 5f 
shell changes. 

II. METHODS 
A. LDA 

We start with all-electron LDA calculation^^^ of the 
band structure of the dioxides at the experimental lat¬ 
tice constants (Tab. [T]). We take into account scalar- 
relativistic effects as well as the spin-orbital coupling. To 
this end, we employ the wien2k packag^i^ with the fol¬ 
lowing parameters: the radii of the muffin-tin spheres are 
i?MT(U) = I?mt(Np) = i?MT(Pu) = 2.65 Ob for actinide 
atoms and i?MT(0) = 1.70 ob for oxygen atoms, the basis- 
set cutoff ATinax IS defined with i?MT(0) x ATmax = 8.5, 
and the Brillouin zone is sampled with 3375 k points (120 
k points in the irreducible wedge). The computed elec¬ 
tronic structure is essentially identical to the nonmagnetic 
LDA results shown in Ref. [HI 

The LDA bands of the actinide 5f and oxygen 2p origin 
are subsequently mapped onto a tight-binding model with 
the aid of the wannierOO codei^ in conjunction with 
the WIEN2wannier interface.!^^ This effective model Htb 
then serves as a base for the LDA-fDMFT calculations. 
The tight-binding representation is not perfect, there 
are some deviations from the original band structure that 
originate mainly in the neglected overlap of the actinide 6d 
states with the oxygen 2p orbitals. The largest differences 
appear in UO 2 , they are explicitly illustrated in Fig. [T] 
It is possible to get a closer match even for UO 2 , but 
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FIG. 1. (Color online) LDA band structure of UO 2 from 
wien2k (thick blue line) and after mapping onto the tight- 
binding model Htb (thin black line). The dominant orbital 
character is oxygen 2p between —8 eV and —4 eV, and uranium 
5f between —1 eV and 2 eV. The visible bands above 2 eV are 
mostly uranium 6d. 

the resulting 5f Wannier functions have a varied spatial 
extension which complicates their interpretation as an 
atomic f shell later on. 


B. LDA+DMFT 

The dynamical mean-field modeling of correlations 
among the 5f electrons amounts to adding a local self¬ 
energy S(z) to the 5f shell of each actinide atom in the 
tight-binding model Htb ■ The selfenergy is taken from an 
auxiliary impurity model that consists of one fully inter¬ 
acting f shell (the impurity) embedded in a self-consistent 
non-interacting medium (i?TB + S) 

The auxiliary model without the f-f interactions can 
be written as 

Hit = E Ploc] Jlh + E Pbath] J\hj 

ij IJ 

+ + ( 1 ) 

iJ 

where the lower-case indices label the f orbitals and run 
from 1 to 14 (or they can be understood as combinations 
of the magnetic quantum number m and the spin projec¬ 
tion a), and the upper-case indices label the orbitals of 
the effective medium that is usually referred to as bath. 
We truncate the bath to contain only 14 orbitals so that 
the local hamiltonian Hioc, the bath hamiltonian Hbath as 
well as the hybridization V are all 14 x 14 matrices. The 
actual determination of Hioc, Hbath and V is discussed in 


detail in appendix]^ Here we just note that they cannot 
be reduced to a diagonal form due to non-commutativity 
of the cubic hybridization with the spin-orbital coupling. 

The truncation of the bath is necessary because the 
Lanczos method, which we employ to solve the impurity 
model, cannot handle much larger systems. In insulating 
oxides, the small bath is well justified on the physical 
grounds: the environment of the 5f shell is dominated 
by the oxygen ligands and hence a small impurity model 
analogous to the ligand-field model should accurately 
represent the dynamics of the 5f shell and its surroundings. 
High accuracy of these small models was demonstrated 
many times in the context of core-level spectroscopies,^^ 
recently for instance in Ref.[TSJ as well as in applications to 
the valence-band electronic structure of transition-metal 
oxides.^ 

The complete interacting impurity model is defined 
as Himp = iitp + U where U is the Coulomb repulsion 
among the f orbitals, 

ijkl ij 

The matrix elements Uijki are decomposed into the ra¬ 
dial expectation values - the Slater integrals Fk, and the 
angular expectation values - the Gaunt coefficients. The 
latter are fully determined by the assumption that the 
angular parts of the local orbitals are spherical harmon¬ 
ics. The Slater integrals Fk are assumed identical for all 
three oxides. They are set to Fq = 6.5 eV, F 2 = 8.1 eV, 
F 4 = 5.4 eV and Fq = 4.0 eV. The average Coulomb 
repulsion U = Fq has a value close to Ref. [TT] to facili¬ 
tate comparison with this earlier study. The remaining 
Slater integrals are chosen to obtain the average exchange 
J = 0.7 eV while keeping the ratios F/^jF^ and Fq/F 2 
equal to their Hartree-Fock values.^ 

The second term in Eq. ([^ is the double-counting cor¬ 
rection that removes the Hartree-like contribution already 
included in the LDA band structure 77x6- In the paramag¬ 
netic phase, the correction can be simplified to an isotropic 
form, = U-a^ij- We approximate Ub by the so- 

called fully localized limit Ub = C/(n/ — 1 / 2 ) —J(n/ — l)/ 2 , 
where rtf is th e self-consistently determined number of 5f 
electrons! ^^ ! ^^ ! The isotropic form of the double counting 
is still only approximate even in the paramagnetic state. 
It is accurate enough for getting the correct position of 
the 5f states with respect to the ligand bands, but it is 
possibly insufficient when it comes to crystal-field effects 
that occur at a much smaller energy scale, especially in 
the case of 4f electrons.!^ 

The selfenergy E(z) in the impurity model i7imp is com¬ 
puted using an in-house exact-diagonalization code that 
combines the implicitly restarted Lanczos method for cal¬ 
culation of the bottom of the many-body spectrurrPH 
with the band Lanczos method for evaluation of the 
one-particle Green’s function.!^ The calculations are per¬ 
formed at room temperature (T = 300 K), well within 
the paramagnetic phase. To lessen the computational de- 
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mands, the Fock space is reduced in a manner analogous 
to the method developed for Ce compounds by Gunnars- 
son and SchonhammerP^ Details of the impurity solver 
are discussed in appendix [B| 


C. Photoemission spectra 


The valence-band photoemission intensity can be evalu¬ 
ated using the Fermi’s golden rule. If the energy and angu¬ 
lar dependence of the dipole matrix elements is neglected, 
the angle-resolved photoelectron spectrum Iy{k,E) is pro¬ 
portional to the one-particle spectral density A{k, E) of 
the tight-binding model, 


A = 


— ImTr 

TT 


E-ir^ 


1 

HTB(k)-f:(E-ir,)_ ■ 


( 3 ) 


A finite imaginary part of the one-particle energy, Fv, 
was introduced to model the life-time broadening of the 
valence states together with the experimental resolution. 
The trace in Eq. (1^ runs over all valence states or over a 
subset of orbitals ^ctinide 5f or oxygen 2p) if an orbital- 
resolved signal is desired. The angle-integrated photoelec¬ 
tron spectrum Iv{E) is proportional to the integral of the 
spectral density A(k,E) over the first Brillouin zone. 

Apart from the valence electronic structure, the impu¬ 
rity model of the dynamical mean-held theory provides a 
means to calculate also the photoemission from core states 
(x-ray photoemission spectra, XPS^^MSo] at the level of 
the so-called multiplet ligand-held theory.^^ To that end, 
we add a core state c and the free-electron states aj, to 
the impurity model. 


.^XPS — T CcC^C-t-C/cv(c^C T ^ , (4) 

k 


where Cc is the energy of the core state, Ck are the energies 
of the free-electron states, fif = fifi is the number 
of f electrons, and Ucv is the strength of the Coulomb 
repulsion between the core electrons and the f electrons. 
The degeneracy of the core state is neglected and hence 
there is only one core-valence Slater integral. 

Employing the Eermi’s golden rule, the probability 
Ic{E) of the emission of an electron from the core le vel to 
a scattering state with an energy E can be written 

/e^lm(0| _ _ ^ „ |0), (5) 

E iFc Eq Cq + Ucv’^f 

where |0) is the many-body ground state of the converged 
DMET impurity model Himp (that is, the initial state 
of the photoemission process), Eq is the corresponding 
ground-state energy, and Tc simulates the life-time broad¬ 
ening of the core state. As in the case of the valence 
spectra, we neglected the energy dependence of the dipole 
matrix elements. In addition, we assumed that the free- 
electron density of states is a constant in the energy 
window of interest. 


TABLE I. Basic characteristics of the investigated dioxides: 
the experimental lattice constant a (Ref. 1321) . the nnmber 
of 5f electrons n/ calculated by a projection onto the Wan- 
nier orbitals, and the computed band gaps compared to the 
experimental values. 



a (A) 

LDA 

DMFT 

gap (eV) 

theory experiment 

U 02 

5.470 

2.8 

2.5 

1.9 

2.1 (Ref.IU 

Np02 

5.432 

3.8 

3.4 

2.5 

2.8 (Ref. mi 

Pu02 

5.396 

4.8 

4.4 

2.5 

2.8 (Ref. mi 


The expression shown in Eq. ([^ applies at zero tempera¬ 
ture, T — OK. The finite-temperature case is analogous to 
the finite-temperature formula for the one-particle Green’s 
function, Eqs. (B3), discussed in appendix [b| Just like 
the Green’s function, the photoemission intensity Ic{E) 
is computed using the Lanczos method. 

In the course of derivation of Eq. (§, the c and 
degrees of freedom were integrated out and hence the 
final-state hamiltonian IJfinai = (.^imp — CAvU/), which 
enters the denominator of Eq. ([^ , involves only the f or¬ 
bitals and the ligand states. Notably, iJfinai has the same 
form as iJimp, only the diagonal elements of Hioc are al¬ 
tered. To make the model more realistic, we scale down the 
hybridization parameters in the final-state hamiltonian 
final, V —?> gV, 0 < g < 1. The hybridization between 
the f shell and the ligand orbitals is reduced due to a 
contracti on of the 5f wave functions when the core hole 
is present.fi^^^ The reduction of the hybridization param¬ 
eters can be incorporated into the hamiltonian of Eq. @ 
by introducing an extra term 

(g-l)(ctc-l)^([V]^,/;&,+ [Vt] jy,) . (6) 

iJ 


This term does not have any transparent physical inter¬ 
pretation but it can still be used to formally check that 
the rescaling V —)■ gV in the final state is compatible with 
the steps that led from Eq. Q to Eq. 


III. RESULTS AND DISCUSSION 
A. Valence-band electronic structure 

The LDA-I-DMFT calculations result in an insulating 
electronic structure for all three investigated dioxides. 
The obtained band gaps are listed in Tab. |T] together with 
the corresponding experimental data. The calculations 
slightly underestimate the gaps, which could be straight¬ 
forwardly fixed by an increased value of the Coulomb U. 
Nevertheless, we get the correct trend in the magnitude 
of the gap across the three oxides even without any such 
fine tuning, simply with the same Coulomb parameters in 
all cases. This indicates that the material-specific screen¬ 
ing of the Coulomb interaction is accurately represented 
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FIG. 2. (Color online) The fc-integrated spectral density corresponding to the 5f states (life-time broadening Fv = 0.4 eV, thick 
red line) is compared with the x-ray photoemission spectra at the aluminum Ka line (connected blue dots). The experimental 
data are adopted from Ref. [33] (UO 2 ), Ref. 1341 (NDO 9 '). and Ref. 1351 fPuO^b The yellow stripes indicate where the band gap 
would appear if the broadening were removed. The Fermi level is placed in the center of the gap. 
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FIG. 3. (Color online) The spectral density of the 5f states in 
the first 2.75 eV below the Fermi level. The life-time broadening 
is smaller than in Fig. (Fv = 0.08 eV) in order to resolve 
finer details of the band structure. 


already in the employed minimal model of the valence elec¬ 
tronic structure. An analogous conclusion was reached in 
a recent LDA-I-DMFT study of transition-metal oxides.l^ 

The fc-integrated spectral density corresponding to the 
actinide 5f states is shown in Fig.j^in comparison with the 
angle-integrated x-ray photoelectron spectra measured 
with incident photons at the aluminum Ka line (1487 eV). 
At this energy, the photoionization cross section for the 
oxygen 2p states is negligible and hence the experimental 
data contain essentially clean 5f signal. The total spectral 
density is measured with incident photons at He II line 
(40.8 eV) as their cross section with actinide 5f and oxygen 
2p states is approximately equal. Comparison of He H 
spectra of PUO 2 with our LDA-I-DMFT calculations is 
shown in Ref. Ell 

Inspecting the case of UO 2 in Fig. [^first, we see three 
distinct features in the occupied part of the theoretical 
spectrum: the main 5f bands near —1 eV, a satellite peak 
near —10 eV, and a broader feature between —3 eV and 
—7 eV. The last feature coincides with the location of the 
oxygen 2p bands and reflects the hybridization of the 5f 
states with the ligands. The three features are discernible 
also in the spectra of Np 02 and PUO 2 , only the 2p bands 
move closer to the main 5f peak, and eventually over¬ 
lap with this peak in PUO 2 . The theoretical calculations 
closely reproduce the shape of the experimental spectra as 
well as their evolution from UO 2 to PUO 2 , which demon¬ 
strates the accuracy of the LDA-I-DMFT modeling of the 
electronic structure. 

The momentum-resolved spectral density is shown as 
a color map in Fig. The dispersion of the 5f bands is 
approximately 0.5 eV in UO 2 and it increases only a little 
to about 0.75 eV in PUO 2 . The hybrid DFT calculations 
predict a much larger increase of the 5f bands dispersion in 
Pu02 due to the degeneracy of the plutonium 5f and oxy¬ 
gen 2p bands (Fig. and a consequent enhancement of 
the effects of hybridization.l^lTlSl Qur calculations suggest 
























5 


U02 

DMFT 

t/=6.5eV 

1 

- 


r* 

Hubbard I 

U=3.0 eV 

1 

¥ 


-10 -5 0 £(eV) 


FIG. 4. (Color online) The spectral density of the 5f states 
calculated with the DMFT selfenergy and 17 = 6.5 eV (red, 
top), and with the selfenergy from the Hubbard I approxi¬ 
mation and U reduced to 3.0 eV (blue, bottom). The only 
larger difference is the absence of the —10 eV satellite in the 
Hubbard I calculations. 

that such increase of the 5f states dispersion is an artifact 
of the single-determinant approximation employed in the 
DFT calculations. Although the angular-resolved photoe¬ 
mission (ARPES) experiments were performed on Pu02, 
the acquired data are not able to resolve the issue yet 
due to the lack of orbital resolution (the measurements 
were done with the He II light source) .1^ The analysis 
of ARPES data is more straightforward in UO 2 where 
the 5f bands are well separated and it is possible to read 
out their dispersion.l^ It comes out as about 0.13 eV, 
which is considerably smaller than our result plotted in 
Fig-i A large part of the discrepancy can be attributed 
to the smearing of the experimental data that merges the 
two bands into a single peak with an apparent dispersion 
approximately 50% smaller than the resolved bands. The 
multiband composition of this single peak shows up in 
the experiment as a variation of the peak’s shape along 
the momentum path. 


B. Screening of the Coulomb interaction 

It was demonstrated in the previous section and also in 
Ref. [TT| that the LDA+DMFT method yields an accurate 
electronic structure if the Coulomb U is set around 6.5 eV, 
or perhaps even larger since our gaps (Tab. |T]) are slightly 
underestimated. On the other hand, the LDA-I-U method 
recovers the correct band gaps in the low-temperature 
ordered phases already with U around 4.0 eV.^ This 
difference originates in the different ways the screening 
by oxygen ligands is accounted for in the two methods. 
In LDA+DMFT, the oxygen states explicitly enter the 
impurity solver as the bath and hence the U used in 
the solver does not include the screening by these states. 
The LDA+U method, on the other hand, is a static limit 


of the LDA+DMFT where a simplified Hartree-Fock 
approximation in the atomic limit plays the role of the 
impurity solver. In this case, the screening by the ligand 
states can enter the calculations only implicitly in the 
form of a reduced U. _ 

The so-called Hubbard I approximatiort^^HU jg analo¬ 
gous to the LDA+U method with respect to the screening 
of the on-site Coulomb interaction. In the same time, it al¬ 
lows for a description of the paramagnetic phase and hence 
can be straightforwardly compared to our LDA+DMFT 
results. The impurity model for the Hubbard I approx¬ 
imation is given by Eq. Q with V = 0 and Hbath = 0. 
Figure shows that the spectral density of UO 2 from 
the LDA+Hubbard I method combined with a reduced 
U = 3.0 eV is almost identical to the spectral density 
from the LDA+DMFT calculations with t/ = 6.5 eV and 
with all other inputs unchanged. It is the screening by the 
ligand states that is responsible for the difference in the 
Coulomb parameter. The same behavior is found also in 
Np02 and PUO 2 (not shown). A related discussion of the 
screening effects in transition metal oxides can be found 

in Ref. [36l 

The close similarity between the LDA+Hubbard I and 
LDA+DMFT spectra of UO 2 indicates that the main 5f 
peak near — 1 eV is a lower Hubbard band. It was argued 
in Ref. |TT] that this feature is a Zhang-Rice state resulting 
from a coupling of the local moment in the 5f shell to 
a hole in the ligand orbitals.l^ But since the Hubbard I 
approximation is clearly able to describe this peak, it 
cannot originate in the Zhang-Rice physics because the 
Hubbard I selfenergy has no knowledge of the ligand 
states. 


C. Core-level spectra 

The nominal electronic configuration of the actinide 
atoms in the investigated dioxides is the 4+ ion, which 
translates to the number of 5f electrons n/ equal two, 
three, and four in UO 2 , Np02, and PUO 2 , respectively. 
Our LDA+DMFT calculations converge to solutions 
with more electrons in the 5f Wannier orbitals, with Uf 
being close to half-integer values, see Tab. |l] A simi¬ 
larly increased filling of the 5f shell was reported pre¬ 
viously in LDA+U , hybrid DFT, and embedded-cluster 
calculations. b^ l ^^ l^The enhancement of the 5f occupation 
indicates a large covalent mixing (hybridization) between 
the actinide 5f states and oxygen 2p states. An evidence of 
such mixing was found also in the experimental core-level 
spectra where the hybridization with ligands is responsible 
for the appearance of satellite peaks.^S An analogous half¬ 
integer occupation of the valence f shell is well established 
in Ce 02 on the basis of the core-level spectroscopy.l^^®^ 
Surprisingly, the earlier LDA+DMFT study of the ac¬ 
tinide dioxides reported an integer occupation of the 5f 
shell in all of them.E^ 

In order to test whether our estimates of ny are com¬ 
patible with the measured core-level spectra, we have 
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FIG. 5. (Color online) Calculated spectra of x-ray photoemission from the 4f core level (red lines) are compared with the 
experimental data (black dots) from Ref. |43] (UO 2 ), R.ef. 1341 ('NdO^I. and Ref. 1441 ('PuOtI. The theoretical and experimental 
spectra are aligned at the 4f7/2 line that is placed at the zero of energy. The life-time broadening Fc was adjusted to match the 
width of the experimental 4f7/2 line in each oxide separately. A background due to the secondarily scattered electrons was added 
to the theoretical curves as described in Ref. 1451 


calculated the photoemission from the 4f states using the 
method outlined in Sec. IIC Since the theory is formu¬ 
lated for a non-degenerate core state, the total 4f spectrum 
is approximated as a sum of two independent components 
(4f5/2 and 4f7/2) that are weighted with the statistical ra¬ 
tio 3 : 4. The splitting between the 4f ^/2 and 4f7 /2 levels is 
taken from the all-electron LDA: 11.0 eV (UO 2 ), 11.8 eV 
(NPO 2 ), and 12.6 eV (PUO 2 ). We employ two empirical 
parameters in the calculations of the core-level spectra: 
the core-valence Coulomb repulsion, Ucv = 6.0 eV, and 
the scaling factor of the hybridization in the final state, 
q = 0.85. We keep these parameters the same for all three 
oxides. Our results are plotted in Fig. [^together with 
the experimental spectra. The satellites are well repro¬ 
duced which indicates that the LDA-I-DMFT fillings Uf 
are indeed reasonable. It is possible to analy ze ind ividual 
contributions to each of the spectral features,l^^^^ but we 
do not enter that level of detail here. 


D. Crystal-field states 


Finally, we discuss the role of the cubic environment 
around the actinide atoms. Its most important impli¬ 
cations are the complex multipolar order in the low- 
temperature phases of UO 2 and Np02,^and the absence 
of magnetism in Pu02.^ The ordered states are outside 
the scope of the present study. The case of PUO 2 was dis¬ 
cussed in our earlier paper - the cubic environment splits 
the 5f shell such that the ground state is non-degenerate, 
and the temperature dependence of the magnetic suscep¬ 
tibility, which appears due to the thermal population of 
excited states, is reduced due to a cancellation of the spin 
and orbital contributions to the susceptibility.^ 

Experimentally, the crystal field can be probed by in¬ 
elastic neutron scattering (INS) that reflects the local 


TABLE II. Splitting of the lowest 5f multiple! by the cubic 
environment. The present calculations are compared to the 
data inferred from the inelastic neutron scattering experiments 
that were performed above the ordering temperature. 




ground state 
multiplicity 

excited levels (meV) 
(multiplicity in brackets) 

U 02 

theory 

3 

161 (1) 164 (3) 

184 (2) 


Ref. [50] 


150 

170 

Np02 

theory 

4 

76 (4) 302 (2) 



Ref. [ST] 


55 


Pu02 

theory 

1 

125 (3) 226 (3) 

319 (2) 


Ref. [52] 


123 

- 


electronic structure at the actinide atoms. Typically, a few 
lowest excited states can be detected. In our theoretical 
description, these states correspond to the bottom of the 
many-body spectrum of the converged DMFT impurity 
model. We compare the computed crystal-field splitting 
of the lowest 5f multiple! with the hndings of the INS 
experiments in Tab. [Ill In fli® paramagnetic state of UO 2 , 
the experiments detect only two excitations.^^ We obtain 
three but the first two of them are nearly degenerate 
and hence they appear as one peak in the experimental 
spectrum. Indeed, this peak is found to split into two at 
the transition into the low-temperature ordered phase.l^ 
In Np02 and PUO 2 , only the first c rystal -field excitation 
is observed in the INS experiments! ^^ * ^^ ! The theoretical 
prediction is in excellent agreement in PUO 2 but it is less 
accurate in the case of Np02. A limited accuracy of the 
present theory is expected due to the neglect of the non- 
spherical components of the double counting correction as 
discussed near the end of Sec. |IIB[ The good performance 
of the theory in PUO 2 is likely to originate in the highly 
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symmetric 5f charge density in the LDA solution^ 

Our calculations yield the same ground-state multiplic¬ 
ity as an alternative first-principles calculation based on 
the LDA-I-U method!^ but the ordering of some of the 
excited states in UO 2 and PUO 2 is different. A direct com¬ 
parison of the crystal-field parameters between Ref. |S3] 
and our theory is not possible since we have non-spherical 
terms not only in the the crystal-field potential inside 
Hioc but also in the hybridization with the ligand states, 
each contributing approximately half of the splitting due 
to cubic environment.!^ 


to the asymptotics of the so-called bath Green’s function 
Go( 2 ;), which is the local Green’s function corresponding 
to the effective medium [Htb + S). The procedure follows 
the steps outlined in Ref. [M] with the notable difference 
that the local hamiltonian Hioc now contains a strong spin- 
orbital coupling that does not commute with the cubic 
hybridization function V(z — Hbath) Therefore, the 

problem cannot be simplified to diagonal matrices. 

The asymptotic expansion of g|q^( 2 :) can be found by 
a repeated application of the identity 

{A-My^ =A-^ (A3) 


IV. CONCLUSIONS ^^ich yields 


We have demonstrated that an implementation of the 
LDA+DMFT method where the selfenergy is obtained by 
the exact diagonalization of a finite impurity model pro¬ 
vides an accurate description of the electronic structure of 
the early actinide dioxides in the paramagnetic phase. The 
band-gap opening is not linked to any long-range order, 
and the main features of the valence-band as well as the 
core-level photoemission spectra are well reproduced. The 
method allows for a quantitative first-principles analysis 
of the covalency between the actinide 5f and oxygen 2p 
states that is found to induce a nearly half-integer filling 
of the actinide 5f shell. 
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Appendix A: Construction of the impurity model 


In this appendix we discuss how the parameters of the 
finite non-interacting impurity model from Eq. 0 are 
found so that the model matches the effective medium 
(the bath) as closely as possible. The hamiltonian has the 
form a block matrix 



/ Hioc 



(AI) 


where all blocks are 14 x 14 square matrices. They are 
determined by comparing the large z asymptotics of the 
local block of the Green’s function,!^ 


^ocW= ^ - Hioc - V(z - Hbath) V 


(A2) 


The asymptotic expansion of the bath Green’s function 
is obtained by expanding the denominator in Eq. (A5), 


Gro{z) - ^2 


yn-\-l 


n—0 


= ye"Ao(e)de, (A7) 


where M„ are moments of the spectral density. The spec¬ 
tral density Ao(e) is a hermitian matrix and hence its 
moments are hermitian matrices as well. The moments 
can be written in terms of contour integrals in the com¬ 
plex plane. Using the path segments depicted in Fig. 
we have 


1 

27ri 




z^Goiz) dz 


1 

27ri 



z"'Go(z) dz , 


(A8) 


that is, an integral over a circle C- r\C+ that encloses the 
entire support of Ao(e). In the calculations reported in 
this paper we encircle the real-axis segment from —20 eV 
to 10 eV. 

When the bath contains the same number of orbitals as 
the impurity shell, the blocks Hioc, Hbath, and V are all 
square matrices of the same size. It is then straightforward 
to construct these matrices by comparing the asymptotic 
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FIG. 6. (Color online) Contours in the complex plane used for 
integration of the moments M„ (shown in blue). Line segments 
are denoted as L±, half circles as C±. The smaller contour 
shown in red is employed to calculate the band occupations 
during the DMFT cycle (the vertical line indicates the Fermi 
level). 


expansions, Eqs. ( |A4[ ) and (A7), term by term. The l/z^ 
terms give the local hamiltonian, 


Hioc = Ml. (A9a) 

The l/z^ terms lead to the equation = M 2 — M^ 
which has a solution 

V = = Y^Ma - M2 (A9b) 

as long as all eigenvalues of M 2 — Mf are non-negative. 
If some of them were negative, the bath Green’s func¬ 
tion would not be representable by an impurity model, 
since the product is always positive as follows from 

> 0 for any jr/j). Finally, the 
1 j terms coincide if we set 

Hbath = V-i[M3-bM3-MiM2-M2Mi](Vl')”^. (A9c) 

Equations ( |A9| ) define all blocks of as hermitian 
matrices. Alternatively, we could take advantage of the 
freedom to arbitrarily choose the basis in the bath segment 
and diagonalize Hbath instead. The transformation reads 
as 


Hbath ^C-iHbathC, V^VC, (AlO) 

where C is the appropriate unitary matrix and the new 
V is no longer hermitian. 


Appendix B: Lanczos impurity solver in a reduced 
many-body basis 

The selfenergy S(z) is computed from the Green’s func¬ 
tion of the interacting impurity model Llimp as 

t{z)=z-ii^^-G-\z). (Bl) 


Since the Coulomb vertex 1] defined in Eq. ([^ only acts 
among the impurity orbitals, the selfenergy is non-zero 
only in the local block Eioc(2). Consequently, it is only 
the local block of the Green’s function that needs to be 
explicitly evaluated, 

Eioc(^) = ^ - Hioc -v(z - Hbath)” Vt - . (B2) 

Simplification to a diagonal representation is not possible 
because Hioc and V(z — Hbath) do not commute. 
The Green’s function reads as 

Gioc( 2:) = X! X! > (B3a) 

N a 

where the two components are 

[G>„(z)]^. = (iVa|/,- ^ ^ . /;|jVa), (B3b) 

Z + h/JSfoi Jlinip 

[^NM., = {Na\f! - \ f]\Na). (B3c) 

The sums represent the grandcanonical average, they run 
over the filling of the impurity model N and over the many- 
body spectrum a, HimplNa) = ENa\Na). The spectrum 
is calculated independently for each filling because a sep¬ 
arate Hilbert space Hn is associated with each N. The 
grandcanonical weights have the form 

WNc = , (B4) 

N a 


where we set the chemical potential to zero without any 
loss of generality. At low temperatures T (large /3 = l/T), 
only a few weights wno. have an appreciable magnitude 
and hence both sums in Eqs. (B3) can be truncated for 
an increased computational efficiency. 

Practically, the bottom of the spectrum including all 
degeneracies is found for each relevant N using the im¬ 
plicitly restarted Lanczos method implemented in the 
ARPACK software package.!^ The matrix elements of the 
Green’s function, Eqs. ( |B3[ ), are e valuated with the aid 
of the band Lanczos method.l^^^^ The band variant is a 
convenient way to access diagonal as well as off-diagonal 
matrix elements of Gioc( 2 ;). 

The Lanczos method for evaluation of matrix elements 
of a resolvent, (z — Himp) , is essentially an expansion 


in powers of Llimp- The length o f the expansion is given 
by the size of the Krylov space.^^SlSfi! The farther is the 
complex energy z from the spectrum of iLimpj the shorter 
the expansion needs to be. During the DMFT iterations, 
we evaluate the selfenergy and the Green’s function along 
the three semicircular contours indicated in Fig. C± to 
determine the non-interacting impurity model and C to 
check the number of electrons below the Fermi level. In 
an insulator, the closest any of these contours gets to the 
spectrum of Himp is half of the band gap (contour C) and 
hence one can get away with a relatively small Krylov 
space. 
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Despite the high efficiency of the Lanczos method, the 
calculations as outlined so far would be very demanding, 
especially on memory, because the many-body Hilbert 
spaces TIn are very large. As it turns out, the Hilbert 
spaces can be substantially reduced without compromising 
the accuracy of the calculations. To this end, we take 
inspiration in the method developed for Ce compounds by 
Gunnarsson and Schonhammer,!^ which can be viewed as 
an expansion in the hybridization parameters V around 
the atomic limit. 

Let us take the ground state of the impurity model in 
UO 2 as an example. The bath orbitals represent the oxy¬ 
gen 2p states and hence they are located several eV below 
the Fermi level (recall the band structure from Fig. [^. 
If there were no hybridization between the bath and the 
f shell, the bath would be completely occupied with 14 elec¬ 
trons and there would be 2 electrons in the f shell. Hence, 
the ground state would be a linear combination of basis 
states 1 /^ 6 °) where b indicates a hole in the bath orbitals. 
If only these states are considered, the bath degrees of 
freedom are completely frozen and the problem is reduced 
to the atomic limit, the so-called Hubbard-I approxima¬ 
tion. When the hybridization is present, the states |/^ 6 °) 
mix with states 1 ^ 2 + 771 ^™^ where m electrons hopped from 
the bath to the f shell. As the hybridization increases, 
basis states with larger and larger m become relevant and 
need to be taken into account. The states can 

be obtained from the zeroth-order approximation \f'^b^) 
by a repeated application of the hybridization part of the 
impurity hamiltonian, 

+ (B5) 

iJ 

hence the connection with the perturbation expansion 
in V. Following the outlined logic, we define truncated 
Hilbert spaces 

= ^\fN-i4+m 5 "*), 0 < m < M} (B 6 ) 

and perform the many-body calculations only in these 
reduced spaces. For the actinide dioxides, we have ex¬ 
plored M up to 4 and have found that M equal 3 is 
already enough to reach essentially converged results for 
all quantities discussed in this paper. The Hilbert-space 
dimensions can thus be reduced from as high as 37 x 10® 
in the case of the full Hn down to about 1.5 x 10® in the 
case of . 

The same truncation of the Hilbert space can be con¬ 
sidered in a general impurity model with bath states 
below the Fermi level and bath states above the Fermi 


level. The zeroth-order Hilbert space corresponding to 
the atomic limit then is 6 ® 6 ®) where b indicates 

electrons in the bath orbitals above the Fermi level. The 
general truncated Hilbert spaces then read as 

Uf'' = 0< m-fn< M}. (B7) 


We used th em in the past to study metallic f-electron 
compounds! ®^ ! ®^ ^ The approximation is variational, the 
truncated Hilbert space turns into the full space 

HN for large enough M. Therefore, although the reduction 
of the many-body basis was formulated as an expansion 
around the atomic limit, it can prove useful also in cases 
where the atomic limit itself is a poor approximation. 


The reduced basis comes with one disadvantage: the 
selfenergy is non-zero also outside the local block Sioc and 
hence Eq. (B2) does not hold. What is worse, Eq. (B2) 
cannot be used even as an approximation because it leads 
to a non-causal selfenergy (alternating sign of the imagi¬ 
nary part of S). Consequently, the selfenergy has to be 
calculated directly from Eq. (Bl| that involves a larger 
Green’s function matrix (28 x 28 in our case). This, in 
turn, necessitates a twice wider band in the band Lanc¬ 
zos method when compared to the calculation in the full 
Hilbert space Hn- Still, the associated increase of com¬ 
putational complexity is significantly outweighed by the 
savings offered by the smaller Hilbert space. 

The appearance of non-vanishing selfenergy outside 
the local block can be under stoo d using the following 
argument. The cutoff in Eq. (B 6 ) can be implemented 
“dynamically” by introducing an artificial multi-body re¬ 
pulsion between holes in the bath orbitals. The extra term 
in the hamiltonian can be schematically written as 


AH 


Upb_^^__^ 

M+1 M+1 


(B 8 ) 


which adds a penalty Up to the energy of states with more 
than M holes in the bath. The strict cutoff is achieved 
in the limit C/p —>■ 00 . The interaction AH apparently 
induces a selfenergy to the bath block, Sbath 7 ^ 0 . 

Recently, a conceptually similar technique for reduction 
of the many-body basis was introduced in Ref. [SH It is 
substantially more sophisticated than the method we 
describe here in that it does not require any a priori 
chosen cutoff M. Nevertheless, it was demonstrated to 
work only in a single-orbital impurity model so far, and 
it remains to be seen how it performs in realistic DMFT 
calculations with multiorbital impurity models. 


kolorenc@fzu.cz 
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